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Abstract. We consider the solution of spectral problems with elliptic coeffi- 
cients in the framework of the Hermite ansatz. We show that the search for 
exactly solvable potentials and their spectral characteristics is reduced to a sys- 
tem of polynomial equations solvable by the Grobner bases method and others. 
New integrable potentials and corresponding solutions of the Sawada-Kotera, 
Kaup-Kupershmidt, Boussinesq equations and others are found. 



1. Introduction 

The paper is devoted to the algorithmic problems associated with integrating 
the spectral problems 

(1) = ■^'i'ix; A) + u^{x, A) -^^^^{x; A) + ■ • • + Un{x, A) A) = 0, 

where Uj are elliptic functions of x and arbitrary (rational or transcendental) func- 
tions of A. We shall restrict our consideration to the Schrodinger equation 

(2) -7^(2:)* = A*, 
the equation 

(3) -u(.t)*' = A^-, 
and the generalization of the Halphen equation 

(4) -u{x)^' -v{x)^ ^ 

We use the term potential for the u{x), ti(x)-fmictions. Until the 1970/80's, few 
exactly solvable potentials were known. Earlier, in 1872, Hermite developed an 
approach for the integration of the Lame equation 

(5) *"-n(n+l)p(a;)* = A^', 
and later, Halphen extended it to the third order equation 

(6) - {n^ - 1) p{x) - i {n^ - 1) p'{x) ^ = \'f. 

Here and below a, (, p, p' denote the standard Weierstrassian functions. See 
[Tsl ^ for an extensive discussion of these classical examples. According to mod- 
ern terminology, the set of exactly solvable elliptic potentials is a particular case 
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of finite-gap potentials in elliptic functions. An intense investigation of elliptic po- 
tentials was initiated by the paper and in 1987 Verdier and Treibich 
unexpectedly found new potentials for the equation (2) in elliptic functions 

(7) u{x) ^ 6p{x) + 2p{x ~ ujj) 

and gave the term elliptic solitons to them. Recently V. Matveev drew attention to 
the fact that such potentials, in Jacobian form, were already considered by Darboux 
in a short note Jill in 1882. The following year, two comprehensive memoires by 
Sparre Jsst appeared on this topic. 

The development of a theory led to the current result that elliptic solitons are the 
widest class of finite-gap explicit solutions. See, for example, |l^, recent results 
in p^ , p6[ and references therein. The monograph ||l| reviews work in finite-gap 
theory up to the beginning of the 1990's and the review |^l[] and preparing book 
|23[ contain a wide bibliography on that score. 

One feature of elliptic solitons is the potential, and the ^'-function can be found 
by the Hermite ansatz method In the case of the potentials with the only 

pole in a parallelogram of periods, the derivation of the algebraic curve and other 
characteristics is not difficult. For this purpose it is enough to take a few resultants 
|l4| 01, but in the general case the elimination technique is insufficient. Sect. 2 
contains a pure algebraic interpretation of Hermite's method. In Sect. 3, we show 
that the general scheme for solving the problem under consideration (1) is reduced 
to the computation of the Grobner basis for some polynomial system. After Buch- 
berger's discovery in 1960's of an algorithm for finding the polynomial ideal bases, 
this area of algorithmic mathematics has rapidly developed. Sec with regard 
to the modern achievements in this area. Sect. 4 contains a relation between the 
method and traditional objects in finite-gap integration theory: the canonical form 
of an algebraic curve A) = and reduction of one of the holomorphic differen- 
tials to the elliptic. Some new examples of elliptic solitons and their applications 
to the integrable partial differential equations (pde's) are presented in Sects. 5-6 
and development of the theory is discussed in Sect. 7. 



2. Algebraic characterization of the Hermite method 
Based on the ^-function 

(8) ci>(x;a) = ^l^eC(")^ 

cr{a) cr{x} 

or more precisely I'element simple 

$(a;;a,fc)^ '^(^ - ■^) (C(o)+fc) . 
a{a) <t{x) 

by Halphen Hermite and Halphen Q considered the following ansatz 

for the ^'-function 

(9) * = $(a;;a,fc) + a^(^' {x\a,k) -\ h a„$^"^ (a;; a, fc). 

The function $(a;;a,fc) as well its a;-derivatives $'(x;a, /c), ... , (x; a. A:) are 
doubly-periodic functions of x of the second kind. According to (1) and (^j-^, the 
expression L^/$(x — XQ]a,k) is a 2-periodic meromorphic function with only one 
simple pole at the point Xq. It must be a constant function. Setting it to be equal 
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to zero, we have L'i! = under the corresponding choice of additional parameters 
fc, a and Uj. As the $(a;; a, fc)-function has the first order pole at a; = p. 231] 

$(.;a,fc) = l + fc+^^^x+^— M4i±^,^ + ... 
X 2 6 

to solve the problem, it is sufficient to equate to zero only the principal part(s) of 

the Laurent's expansion(s) of the expression L^f, where 5* is the ansatz or its 

multi-pole generalizations (see Examples 1—3, 7 in Sect. 5). 

As a nontrivial example we shall consider the 5-gap Lame potential with n = 5 
in (H). It has been studied in |Q, but we give a more simple solution. Note that 
the cases n = 2, 3, 4 and partially 5 were considered already in ||2^, pp. 527-531]. In 
the same place one can find mention of eliminations. 

By Frobenius theory, ^ must have a 5-th order pole at x = and therefore the 
ansatz for the ^f-function should be the following: 

(10) * = $ + a^$' + ---+a4$(''''. 

Substituting (HO) in (H) and expanding the result at a; = 0, we obtain a system of 
equations in the variables k, a, a^. This system is linear with respect to aj. We 
do not write expressions for the Uj (Jli], formula (3.7)]). The remaining equations 
have the form 



xui = -6 fc5 + Y (9 P + A) fc» - 60 p' fc2 + (^90 p2 - 20 A p - y + — j k - 
-|(9p-5A) p', 

W2 = -5k^ + {75p + 5 A) - 100 p' + (225 p^ - 30 A p - ^ + y^ 32) fc^ - 

- 20 (3 p - A) p' fc -f 25 p^ - 15 A p2 + I (A^ - 20 52) p - 40 (73 - ^ A=^ + y ff2 A, 

which are understood to be equal to zero. The argument a in the p, p'-functions is 
omitted for brevity. The system (|^) has to be considered as algebraic with respect 
to k and transcendental in a. We emphasize that everywhere in the paper A is a 
parameter, but not variable in polynomial bases. Insomuch as functions p(a) and 
p'(a) are related by the Weierstrass equation (torus) 

(11) Wr^ = p'{a)^ -Apiaf +g^p{a)+g^, 

we supplement (^ by ( pi] ) and consider (|^ 11) as a polynomial system with respect 
to independent variables {k, p, p'). The simplest method of solution consists of the 
elimination of the variable p followed by p'. As a result, we find that k must be a 
root of the polynomial 

k^ (5103 k'^ - 945 A fc^ + 40 A^ + 54 g^f (225 (27.92 - A^) PeW + -Pgl-^) -P|(A)), 

where Pg g 3(A) are some polynomials in A with degrees 6, 9, 3 respectively. It is 
not difficult to guess that the correct result requires that k and A are related by 
the following equation 

(12) F(fc,A): 22h{27g^-\^)Pl{\)e + P,{\)Pl{\)=0, 

since the differential equation (|^) is of second order and we must have not more 
than 2 solutions for k with fixed A. Curve ( |l2|) can be brought into the canonical 
hyperelliptic form 

F(/x,A): /i2 = (27g2-A2)P9(A) 



4 



YURII V.BREZHNEV 



by an obvious birational transformation (note a misprint 27 in this formula in 



14, formula (3.8)]). The variables p, p' as functions of k can be found in the same 
way: by the sequential reduction of exponents of p, p' in d^-pl]) . 

Obviously, the resultant technique is almost impossible if the potential has sev- 
eral poles |3^ , as the number of variables increases. Another approach consists 
of finding an equivalent system with the following criterion. It is advisable for the 
new system to contain linear equations in p, p'. These equations define a as a 
function of (fc, A): 

(13) p(a) = i?i(fc; A), p'(a) = R^(k; A) 



and we call (13) a cover of torus ( |l l| ) in algebraic form. Suppose one of the new 
equations does not contain p, p' (i.e. be an univariate polynomial in fc) if the 
nontrivial solution for k exists. We interpret such a polynomial as the algebraic 
curve F{k, A) = corresponding to an elliptic potential. If F{k, A) has a factorised 
form then each of factors is investigated separately. The curve is one of them. 
It is clear that its degree in k has to be equal to the order n of the equation 
(1). The canonical form F(/i, A) = of the curve is obtained with the help of a 
birational transformation between variables (k, A) <-> (/x, A) (see an explanation in 
Sect. 4). Note, there are specialized algorithms for the computation of the univariate 
polynomial in an ideal without solving the system as a whole. 



3. Grobner bases, curves and covers 

We clarify the main idea using the previous example. Let us consider three 
polynomials Wi 2 3(p' P: fc) as a system generating an ideal in a polynomial ring 
Q(A,52>53)[p'>P> k] 

(I) ^ h^w^ +/12 W2 + 
where — hj{p', p, fc) are arbitrary elements of the ring. As is well known, the 
structure of the solution of the polynomial systems depends on the monomial or- 
dering in a ring pl| |. The arguments at the end of Sect. 2 (see also the elimination 
theorem in |10|) lead to the choice of pure lexicographic ordering p' ;^ p ;^ fc. The 
monograph [|lO| contains a good exposition of details of this subject. The main 
property of the Grobner base is expressed in 

Definition Let {wi, W2, ■ ■ ■} be a basis of ideal I = {wi, W2, . . . ). Let >~ 

be a monomial ordering on the ring Q[. ..] and LT(/) denote the leading term 
(monomial) of a polynomial f e Q[. . .]. The set G = {fi, /2, . . . , /jy} is said to 
he a standard basis ( Grobner Basis ) if the monomial ideal generated by 

(LT(/i),LT(/2), ... ,LT(/^)) 

is coincident with an ideal (LT(/)) generated by all the leading terms of L. 
In other words, the leading term of any polynomial in / is divisible by one of the 
LT(/^). According to the definition at the end of Sect. 2, the polynomial F{k;X), 
determining the algebraic curve 

(14) F(fc;A)=0, 

is a generator of the intersection of the ideal / and the ring of all polynomials in k: 

{F{k;X))=inQW[k]- 
Thus we arrive at a general recipe for the solution of the spectral problem (1). 
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Proposition. Let {wi, W2, ■ ■ ■} be polynomials in p'(a), p{a), k, . . . appearing in 
the Hermite method and determining the solution of the spectral problem (1) : the 
cover of torus ( p^ ) and the curve (|l|). Then 

1) Algorithmically, the method of solution is reduced to the computation of a 
standard basis for the ideal I = {wi, W2, • • . } with respect to pure lexicographic 
ordering p'{a) >- p{a) >- k >- .... {For example, by Buchberger's syzygy 
polynomials algorithm Q); 

2) Let G = {/i, /2, . . . , / jvj be this basis. The algebraic curve (Q) and its 
projection on the torus ( |ll[ ) in algebraic for m ([13| ) are contained in G if the 
univariate polynomial in k and polynomials ( |l3|) exist; 

3) IfG contains a polynomial free from variables fc, p'(a), p{a), then the spectral 
problem (1) is not integrable in the framework of Hermite 's ansatz. 

Proof . 

1) The standard basis always exists and Bucliberger's algoritlim terminates pO| . 

2) Taking the resultants oiwi, W2, ■ ■ ■ we eliminate variables p, p' and get poly- 
nomials) R{k). It is obvious that R{k) g /. Using the divisibility 

(LT(A),LT(/2), ... ,LT(/^)) = (LT(/)) 

and lexicographic ordering, the equality R(k) = h{k) f has to occur for some f ^ G 
and h{k) S Q(A)[fc] (possibly equal to 1). Therefore there exists a polynomial 
/ depending only on k. Designating f{k) = F{k;X), we obtain the curve (|l4|). 
If F{k; A) has a factorized form, then the algebraic curve is one of its factors. 
Analogously, if the polynomials (|l3| ) exist, then they necessarily belong to G. In 
the same way, an important formula — the curve as a cover of the torus (|l^) in a 
transcendental form (an equation in a) 

(15) i?(p'(a), p(a); A) = 

necessarily must be contained in G computed with the ordering k >- (p'(a) >- 
p{a) >- ■ ■ ■), where permutations inside the brackets is allowed. Note, the order of 
elliptic function (|l5| ) in a is equal to the order n of the equation (1). 

3) An existence of such a polynomial implies a restriction on the spectral pa- 
rameter A (see a demonstrative Example 5). ■ 

Note a direct link of the point 3) to a treatment of finite-gap potentials as Picard's 



potentials |20 



There are numerous algorithmic methods to solve this problem. Among them: 
the Grobner basis method ||, the method of characteristic sets, and an effective 
method of elimination based on the Seidenberg theory We do not discuss all 
the modern achievements in this area. See and references therein for details. 
Note that the reduction of the holomorphic differential da to the elliptic one is 
derived from ( |l3| ) by the formula 

^ ' P'(«) R2F, 

where subscripts fc, A denote the derivatives with respect to k and A. 

4. Canonical form of curves and holomorphic differential 



Formulas (14-16|) give a noncanonical form of the curve and holomorphic dif- 



ferential, i.e. expressions in the variables (fc. A). The canonical variables (/i. A) in 
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the algebraic curve F{fi, A) = we call variables A in (1) and eigenvalue /it of a 
commuting operator pencil 

(17) P(A)^' = M^'- 

Supplementing the polar expansion of the equation (1) by the polar expansion of 
(P^, we get the algebraic equations in variables (p, p' , fc, /i, . . . ). Again, based on 
the above properties of the Grobner base, a canonical representation of the solution 
and all the spectral characteristics are extracted by computing the base with the 
ordering (p >- p' >-...)>- fc ;^ /i. Such a base contains a birational transformation 
between the (fc, /i)-variables in one direction: 

(18) fi^k: fc = i?3(/i;A). 
An inverse transformation 

(19) fc^/i: fi^R^{k;X), 

where 4 are rational functions of its arguments, is computed by the ordering 

(p ^ p' ^ . . . ) ^ ^ fc- 

We add a few words about the efficiency of computations. The solution of a 
spectral problem itself does not require the inclusion of a commuting operator 
(p^l). So, among of its polar expansions one may take (and supplement) those, 
containing only the /i-variable. Evidently, it will enter into the polar expansion 
with first degree: 

(20) ii = wik,p,p';X). 



After the computation of the base (not including (O)), we will have the curve (14) 



and cover in algebraic form (y_3|). Substituting it into the equatio n (|20| ), the pair of 



equations for the determination of the above transformation (|lg, |19| ) is 

(21) ^ = u;(fc, i?i(fc; A), i?2(fc;A); A), F(fc; A) = 0. 

Formulae (|l^, |l^) are obtained by computation of the bases for (^ij) with ordering 
(fc >~ /i) and (/i >- fc) respectively. See Example 6 for details. 



5. Examples and applications 

In this section we demonstrate the ideology of Sects. 3-4 on examples. The 
generality of the technique allows us to make further proofs. Let us prove that the 
well known potential of Treibich and Verdier (0) |Q for the equation (2) is the 
only possible 2-pole potential in the class 

(22) u{x)^6p{x) + 2p{x-^}), n^O. 

Example 1. The Treibich Verdier potential. Parameters A, g2, 53 are 
fixed and f2 is an unknown constant. The ansatz for the ^'-function must be the 
following: 

(23) * = aQ^{x; a, k) + a^^^' {x; a, fc) + a2<i>(a; - Vl\ a, fc). 
Substituting (|2^, |2^) into (2) and equating the poles to zero, we obtain ^'-function 

5- = 6 $(-rj; a, fc) $'(a;; a, fc) - (3 fc^ - 3 p^ - 2 p^ - A) $(x - VL- a, fc) 
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and a system of five polynomials: 

= 2(p„-pjfc3 + 3(p; -p^)fc2 + 2(p„-pJ(3p„-A-2pJfc + 

+ (6 - Pa) Pa - (Pn " Pa) " ( ^ Pa " 2 P^ ) P^ 

= 3/c3-(9p„-4p^+A)/c + 3p^+3p;, 

(24) w, = 3/c4-2(9p„-14p^^-A)fc2 + 12(p/^ + p:)fc- 9p2- 

- 2 (14 p^^ + A) p„ + 12 p.^' - 8 A p^^ - A^ 

"^4 = pf - 4 p_i^ + 52 Pa + .93 

^5 = pf -4p„^+52Po+53- 
where we used the addition theorems for elliptic functions, the important equality 

$(17; a, k) $(-J7; a, fc) = p„ - p^ , 

and designated p^ = p(a), p^ = p(ri), etc. A common factor p^ — p^ was removed 
in polynomial because it leads to the contradiction: a = (aj — oo). The 
system (|4|) generates the ideal 

(25) W2, W3, W4, W5) G Q(A,g2,53)['i;>Pa'Pa>Ps^>Pn]- 

Note, the same ideal (wi, . . . , w^) in the ring Q(A, g2, g^, pj[k, p^, p„, p^] leads to 
the just mentioned condition p^ — = 0. Therefore is not arbitrary. Computing 
the minimal reduced Grobner basis for (|2^) with pure lexicographic ordering p^ >- 
Pa ^ ^ pQ ^ Pq: obtain 8 polynomials, some of them having a factorized 
form. If some of the factors do not depend on A, we obtain restrictions on p^ and 
p^ equating these factors to zero. There are four such polynomials: 

Gi = p^ ((A3 -4.92 A- 16 <?3)/c+ (3 A2- 4^2) P^)M, 

G2 = p^ (l6p;A; + 3A2-452) M, 

G3 = (4p^^' -52p, -53)(4p„ -A)M, 
G4 = p^(4p„-A)Af, 

where the multiplier M denotes 3 A:^ — A — 5 p^ . The equation M = yields the 
trivial result p{a) = p{^)- It is checked by recomputing the base ( p^ ) with an 
additional polynomial M. Further, fl must not depend on A(!). Therefore, the 
only solution for Q is defined by the equation 

p'{n) = Q =^ n = ujj^, 
where lu^ are the half-periods of elliptic functions. Substituting p^ = 0, p.^ = Ci, 
52 — 4 (e^ + 65^62 + 62), .93 = — 4e]^e2(e]^ + 62) into (24) and recomputing the basis 
with respect to the ordering Pa Pa >~ k, we obtain the well known algebraic curve 
of genus 2 and all algebraic-geometric objects For classification results of 

the Treibich-Verdier potentials and other elliptic ones, see [^ , the appendix in 
(Tsl and the most recent results in the review 

Example 2. As a preliminary, we shall consider equation (3) with potential 

(26) u{x) = 6 p{x) + 6 p{x - fl) 

and the restriction 52 = 0. As before, we have the ansatz for the ^'-function 
* = $(x; a, k) + a^<^{x ~ VL; a, k) 
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and the original basis of the ideal is generated by 5 polynomials. Computing the 
Grobner basis G, we obtain a system of 8 polynomials. Only two of them have a 
factorized form. 

Gi - (64A3fc3_27(A2 + 16.g3)2)(|:,;fc-3p,f), 

G2 = (64 A^fc^ _ 27 (A2 + 16 g^^f) ((4 - .93)^ - 3 p^') . 

The nontrivial solution will take place if and only if (fc, A) be coordinates of the 
algebraic curve which is the first factor in Gi, G2. In the next example we rule out 
the condition = 0. 

Example 3. If ^2 is free, the straightforward computing of the basis is un- 
successful. Indeed, the Grobner base method is universal and therefore it can be 
ineffective in some special cases. But our interest is only with the zero structure 
of the polynomial system. Thus, the characteristic sets method ji^ is the best ap- 
proach in this case. Under the ordering p^j -< pj^ ^ ^ ^ Pq, ^ p^, the characteristic 
set has the form 

A = (64(A;A + .g2)(fcA-252)'-27(A2 + 1653)')A:M, 

/2 = (8(A;A-2 52)p„-8A;3A + 16.g2A;2 + 3A2+48 53)- 

(4p;^fc + 52-12p„')fc, 
/3 = (32(fcA-2 52)A^P:-192 5|fc'A2 + 
(27) + (A* - 288 53 A^ - 6912 gj + 2^g^)kX- 

- 52 (11 A^ + 864^3 A^ + 6912 gl - 2^ gl)) k M, 



where 



M = 64 p; p„ fc2 _ 4 (3 p;^ A - 16 52 P, - 12 53 + 96 pf ) k 



+ 3(12p,^2-52)(A + 4p^). 

Factorisation shows that the variable f2 is separated in polynomials /i 2 3- Therefore 
fi gives an algebraic curve independent of Q: 

(28) 64 (A: A + 52) A - 2 52)^ 27 (A^ + 16 53)^ 

The polynomials {/2, f^} are an algebraic form of the cover (|l3|). However, the 
genus of the curve (|8|) is unity and we have a cover of a torus by a torus. Hence, if 
moduli of both tori are equal, then there is a one-to-one correspondence between the 
global parameter a of the torus and the global parameter r of the torus ( p8| ) . 
The next step is to find it. After the birational change of variables {k, A) <-> (y,x): 

3 -I- 2 On 2; -I- 3 03 , 

k=— X = Ay 

Ay X 

we obtain the canonical form of the curve ( p8| ) as = 4 a;^ — 52 ~ 53 with an 
obvious uniformisation and the equality a = 2 t. The final solution of the problem 
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(3, ^6|) is as follows: 

^-(0;; A) = a 2 T, fc) + $(^2; 2 r, fc) $(a; - 17; 2 r, fc), A = -4 p'(t), 

^^^^ a = C(2r + n) -2C(t) -C(f7), A: = 2 C(t) - C(2 r). 

The passage to the limit r ^ a;^ in ( p9|) leads to the solution under the condition 
A = 0: 

^-(2;; A = 0) = Ci (((2^) - ((^^ - ^)) + C2. 

An attempt to integrate the more general potential 

u{x) — 6 p{x — ill) + 6 p(a; — ^2) + ^ 

with a nonzero constant A failed. However, this point has an explanation in the 
theory of nonlinear partial differential equations. Indeed, the spectral problem (3) 
is associated with the Sawada-Kotera (SK) equation 

(30) u^ — ^3,3,3,3,3, — 5u 7/3,3,3, ~ 5 M3, u^^ + 5u'^ u^, 

By assuming that the poles Qi 2 depend on time t, one obtains an isospectral 
deformation of this potential. This simple calculation yields the stationary solution 
of® 

u{x, t) — 6 p{x — ct) + 6 p{x ~ ct — Q) + A 

with the conditions 

c+12.g2 + 5A2 + 60Ap(f7) = 0, Ap'{n)=0. 

Therefore (Tl = cuj and A is free) or (A = and is free) . In the both cases we 
obtain a restriction on a velocity c of two cnoidal waves. See an example in [ ^9[ for 
the case A — 0. Recently, Conte and Musette obtained a similar result ||3^, formula 
(84)] and revealed a remarkable more general solution in an old paper of Chazy 
in the context of the Painleve analysis: 

(31) u{x,t) = 6p{x - ct-n; g2,g3) + 6p{x~ ct-il; 02,93), 

c = -1552' c = 352 -1552- 
Strictly speaking, Chazy's solution ||^, p. 380] corresponds to the stationary equation 
( ^0| ) and therefore to the case ^2 = .92 = c) in (|3l|). One can show that the 
potential (^ij) is the stationary solution of a linear combination of the equation 
50|) and higher SK-equation of the 7-th order 



(32) . _ ^ _ _ 3 4 3 \ 

T f ^53; ^ 2 ^4^x 3 ^xx'^xxx 2 ^ ^xxx ^ ^ ^a; '^xx ^x g ^ ^x J ' 

Sect. 7 contains additional information for this potential. We do not enumerate 
other 1-pole elliptic potentials u{x) = A p{x) + B for the equation (3). For example, 
one of them is m = 30 p{x) ± 3 \/3g^ (see also |H). 

Example 4. Let us consider a general 1-pole elliptic potential for the equation 

(i 

(33) 'f'" - (a p{x) +d)<f - {b p'{x) + c p[x)) = A * 
in the framework of the ansatz 

= cf>(a;; a, k). 

Using the above techniques in the ring Q(a, 6, c, d)[p^, Pqj ^] we do not get the solu- 
tion: I — (1). Therefore (a, 6, c, d) have to depend on each other. After calculations 
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in the ring Q(a, 6, c)[p^, p^, fc, d] we determine step by step the constants (a, b, c, d) 
and get the foUowing. The first polynomial in a base is 

-8 fo^ (2 6 - 3)3 - 4 c (2 6 - 3)2 (4 d &3 _ 12 fo2 ^ _ 6 + 6 c^) A + • • ■ = 0. 

Equating to zero the coefRcients in front of A^, A we obtain 

b^-, or 6 = 0. 
2 

In these cases we will have respectively 

Sd-c^^O, (216A + c3 -36cd)c3 = 0. 
Therefore {b = 3/2, d = 0^/3) or 6 = c = 0. In the first case we have 

- 3 {p{x) + c^) - (I p'{x) + 3 c p(x)) = A A) = a, c). 

The nonramified cover of the torus (^5|) of genus g = 1 is 

p'(a) -6c p(q;) + 2 a + 4 = 0. 

c is an arbitrary constant and the condition g2 — ^ [ p8| , example 3.10] does not 
appear. Note, there is no such restriction in the Halphen equation (o) with n = 5 
as in jlj]. It appears only for n = 4 jl^. The second case is known 12^ : 

(34) - (6 p(a;) + rf) = A 

108Ap'(a) + 36(^2 -3 52)p(a) + 27A2 - lOSgg -4^^ = o (genus g = 2). 

See 1^ for an application of this potential. 

Note that both cases can be found in t. Ill: pp. 372, 522] in Jacobian functions 
and jl8[ III/IV: pp. 460, 462] in Weierstrassian functions. No other possibilities ex- 
ist. The same technique is applicable to other ansatzs. The next one is a nontrivial 
example along these lines. 

Example 5. The equation (p3) in the framework of the ansatz 



(35) * = ao$(a;;a,fc) +$'(a::;a,fc). 

As a consequence of corresponding indicial equation ( ^7| ) with = 2, without loss 
of generality we get b = 12 — a in (|33|). Solutions for a, c, d must not depend on 
A and k. One solution suggests itself. Indeed, the first polynomial in the original 
base has the form 

(a - 12) ((a - 18) (fc^ _ p^) + 2 c fc) + 2 (a - 18) d + = 0. 

With a — 12, this polynomial does not depend on fc, a and we get (after the 
replacement c ^ 12 c) 

a ==12, 6 = 0, d=12c^ ao = fc-2c. 

Moreover, the ideal in the ring Q(A, c, 53) [Pq,, p^, k] is not equal to (1) and there- 
fore, c is an arbitrary constant. Thus, the equation (|3^) and its solution take the 
form 

- 12 (p(x) + c2) - 12 c p(x) * = A 

(36) 

^{x; A) = $'(a;; a,k) -2 c^{x; a, k) . 
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We do not give here the large formulae for the cover ([l3|), or the 4-sheet cover in 
the form (^5|) and write only a skeleton of the non-hyperelliptic trigonal algebraic 
curve ( p^ ) of genus g = 3 

64 (A^ + 32 A + 2^ c*^ - 108 53 c^)(A - 11 c^) fc^ + (■••) fc^ + (•••) fc +(••■) = 0, 

where (• ■ ■ ) designate some polynomials in A,c, 53 with integer coefficients [Q. 
Under c = we arrive at the case ( p6| ) with = 0. 

The higher ansatzs (||) are investigated in a similar manner. Indeed, by Frobenius 
theory, if 5" has the expansion = x~'^ + ■ • • , then a, b satisfy the determining 
equation 

(37) + + 2) + aiy + 2b = 0. 



A natural question appears: under what parameters (a, b) does the equation (37) 
have integral solutions for One of solutions is Halphen's equation (||). It corre- 
sponds to 2 6 = a and ( ^T\ ) is reduced to 

a — ~ 1 {n = v + 1). 

As in the previous example we can list all the possible cases for the ansatz (^5|). 
Indeed, assuming & = 12 — a and (a, d, c) to be arbitrary, the origin base contains 
three polynomials The first and second of them are linear in p, p'. Solving 
them and substituting into the base again, we obtain the remaining polynomial in 
(A, fc): 

-4 (a - 12)2 _ -^g)2 ((a _ 6) (a - 18) fc + c (a - 9)) A +(•••) fc +(•••)== 0, 

where dots denote a polynomial in (a, c, d, §2, g^)- It must be zero for all values 
of A. Splitting it in A we get two linear polynomials in k. Their compatibility 
condition is the polynomial 

(a - 6) (a - 8) (IO8 c"* - 72 (a - 18)^ dc^ + {a- 18)'' (12 d^ - {a - 12)^ g^)) = 

and solution for k 

k^ (a-9)c 



(a - 6) (a - 18)' 

The verifying of Weierstrass's relation jll] ) yields a polynomial in A 
(a - 6)3 (a - 12)3 _ -^g^e + {...) \ +{...) = Q. 

Under a 7^ 6, 12, 18 we arrive at the point 3) of the Proposition. Therefore, only 
three possibilities exist: a = 6, 12, 18. The corresponding final solutions for the 
variables (p', p, k) are obtained separately: by recomputing the base. Thus, besides 
(p6[), we have the following integrable potentials (note a misprint p{x) instead of 
p'{x) in one of the formulae in j^) 

^'"'-(18p(a;) + <i)*' + 6p'(x)1' = A*, - {6 p{x) + d) ^' -6 p'{x) ^ = A*. 

Sec for solutions of the generalized Halphen equation (||) and Q for details of 
the Example 5. 

It should perhaps be noted here that the example (^) is the generalization 
c 7^ of the first nontrivial case n = —3 in a series of other Halphen's equations 
lip. 554] 

(38) w'" -\n^w p{z) - ^ n (n 3) (4 n - 3) w p'(z) = 
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without a spectral parameter^]. Notation as in |T8[ III/IV: Ex. 15, p. 464]. Indicial 
equation (B^) for the example (pSh becomes 



(3j/ + 2n)(3i/ + 2n + 6)(3i/-47i + 3) = 

and (n + 3)(n + 6)(4n - 9) = for the ansatz (||) (j^ = 2). 

Example 6. Halphen's equation (||) with n = 5. Here we display only the final 
formulae in the context of Sect. 4: 

• The commuting operator pencil (pT|): 

A VI/" - 14 (4 pixf - 52) + 16 (7 p'{x) - A) p(x) * = ^ 

• The polynomial @: 

6At-56fc^ + 560p„fc^ - 20(28 - A) fc'' + 168(5 -32) - 4 (28 p^ + 5 A) p„ = 0; 

• The birational transformation ( p^ , p^ ), which is quadratic in (fc, /i): 

32 (2 ( A2 - 392 .93) fc A + 7 (5 A^ - 784 33) g^) ((A^ - 392 g,)k-2lg^ A) 
fi — — 



49 A4 - 208 53A2 + 3136 (5I + 4 5^) 

7 m'-452 (5 A2- 78453)^ 



2 

/ / ■ 

k ^ 



(A2- 392 53)A 

• The canonical form of the algebraic curve of genus 4 (see also p8[|): 

F(Ai, A) : _ ^2 _ 784^^) ^ _ + 20833 A' - 3136 {gl + 4gl) A = 0; 

• The 8-sheet cover in the form (p^: 

2** (A^ - 392^3)3 A p'(a) - 2^49 32 (A^ + 112.93) (A^ - 392.93)^ p(a) + A** = 0. 



Note that both this cover and its algebraic form (^_3|) are the expansive expressions, 
whereas the reduced holomorphic differential (|l^) in the variables (A;, A) and (/x, A) 
is given by the simple formulae: 

dpja) ^ -8(A^-56ff3) ^ ~7(A^-56g3) 

p'{a) 3^2 -452 (11 A2- 78453) (A^ - 39253) (3 Afc + Wg^) ' 

Analogs of the above formulae are derived for all other examples in the paper. 
Example 7. The 2-pole potential for the equation (jj) with ansatz 

\E' = <i>(x; a, k) + <i>(x - il; a, k). 

The general 2-pole elliptic potentials contain many parameters — multipliers before 
the p', p, ^-functions. We do not give their exhaustive classification and consider 
only the most interesting case 

- 3 (p(x) -I- p(x - r!) - p{A)) - 

(39) o / ^ 

- { + p'ix-n) + B pix) - B p{x -^))m^X^. 

By virtue of the Proposition, the Grobner base contains all the information about 
the solution, i.e. all the following formulae. As before we obtain that f2, A are 
arbitrary constants and 

B = 2^p{n)-p{A). 



^The example (BSl) was revealed by E. Previato 
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The parameters k and A as meromorphic functions are related by a algebraic equa- 
tion of the genus 2 independent of fl (compare with 

(40) 2Xk'' + i3p^-g^)e-3p^Xk^\x^ + pf = 0. 

The equation (^o|) can be realized as a 2-sheet cover of a torus in the form (^5|) 

A p'ia) + (3 p^ - 52) P(a) + \>^^ + Pa~93^ 0. 
The algebraic form of cover (nSf) has the form 



p{a) ^e + p,, p'{a) = ~ 3 p^) y - X ~ r 

6. Solutions of integrable PDE's 
The spectral problem (^, |3^) corresponds to the Boussinesq equation 

(41) 3u,t^ {2u^~u^^)^^, 

and the arbitrariness of SI means the existence of an isospectral deformation of the 
potential 

(42) u{x,t) ^ 3 p{x - n^it)) + 3 p{x - n^it)) ~3p^. 

Substituting the 2-gap ansatz ( ^ ) in (^Tj) , we get the well known system of pairwise- 
interacting particles of the Calogero-Moser system type with a repulsion 

potential and immovable center of mass. Integration leads to an equation for Sli(i): 

Using the uniformisation of the corresponding elliptic curve, we obtain the explicit 
solution 

(43) ^,{t) = \p-' (C(r + v)-aT-v)-a2v) + \p,) + '-, 

- = ±^r^(| + |), p'(2.) = -^(pi-52a+253), 

where r = 8 (t — tg) and c, tg are arbitrary constants. An implicit form of this 
solution with c = in terms of Jacobi's sn, cn, dn-functions is given in and 
in the earlier citation |^ in the context of solutions of the Kadomtsev-Petviashvili 
equation. In (43), the elliptic integral p~^ is calculated with invariants 53, and 
the C,, p~^-functions with invariants 

2 , Q„ , , 9l\ ~ _ 53 / /2 , \ , .92 52 52 



52 = 77 52 P. 



o , , 52 \ ~ _ 53 / /2 . r, 

"-A ' •^.93 Pa + I : 53 - ypA + ^53y ' 2733 24 



The reduction case _B = in ( p9D corresponds to the Kaup-Kupershmidt (KK) 
equation 

(44) Uf. = u^^^^^ — 5u u^^^ — i^^x ^xx + 5 Uj. 
and its stationary solution 

(45) u{x, t) ^3 p{x - c<) + 3 p{x -ct-n)-3 p{n), 
but the velocity c depends on the distance between poles: 

c = 352 -45 p^{n). 
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Pi - P2 
Pi - P2 



The generalization of ( psj ) in a similar way as the solution ( pl| ) is 

(46) u(a;,i) = -12(pi + p2) + 3 
where 

EE p(a;-ct-f]; g2,53)> p2 = p(a^ - - fi; .92- Sa). c=-12g2- 

Here f2, fi, 53j 53 ^'^^ arbitrary constants. Using a connection between the 
SK(u)- and KK(t(;)-equations and the Tzitzeica equation 

(47) </,,, = e^-e-2^ 
with the Fordy-Gibbons equation 

via the Miura transformations 

we obtain stationary solutions for the w-function 

pi - P2 

Non-stationary solution of the equation ( p7| ) has the form 

(j){x,t) = \n2 c + \n[p{x + cA -n\g^,g^) ~ p{x ~ ct - h; g^,^^)) 



with the restriction: 4(^3 — g^)c^ — 1. The details of calculations (43) are ex- 
pounded in Q and the formula (^), ^'-function for the potentials ( |31| , |4^ ) in 

7. Concluding remarks and discussion 

The investigation of elliptic solitons can be automated by a polynomial tech- 
niques. The Grobner base method provides an unified approach to the solution of 
related problems. The technique suggested with minor modifications is extended 
to matrix spectral problems. 

As the Examples 2—3 and 5 show, the algebraic curves can be degenerate 
(multiply roots of discriminant). 

The general case in Example 7 for the equation (|^) 

u{x) = a p{x) + bp(x — Q) + c({x) — c(^{x — 57), 

(48) 

v{x) — dp'{x) + e p'(x — il) 4- / p(x) + g p{x — Q) + h ({x) — h ({x — V.) 

requires additional research. 

To all appearances, the exampl e (|3^) h as to fit into the hierarchy of higher Boussi- 



nesq equations, studied in full in ^2, 13 . Multi-pole potentials are investigated by 
involving the addition theorem for the ^-function: 

, , s 1 ^{x; a, k) <S>iz] a, k) fp^ + p^ + p'^ 
9[x + z]a,k) — - ' 



2 Px-Pz \Pa-Px Pa- Pz 

A natural assumption suggests itself: all the potentials, obtained by the above 
method, are finite- gap ones. At least, by construction, all such potentials belong 
to the set of exact integrable (explicit ^ |39j) and the ^P-function is a single- value 
function on a Riemann surface of the algebraic curve F{k, A) = and meromorphic 
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function (in x) for all values of A (Picard's theorem ^). Note that remaining 
linear independent solutions for the ^'-function are got by choosing of fc-branch of 
algebraic equation (p^). 

If the assumption is valid, then the potentials for the spectral problems (2-^ are 
free of residues (a consequence of 0-formulas), and therefore ansatzs for the multi- 
pole potentials (say (^)) do not have to involve the functions. This strongly 
decreases the number of parameters and the computational task. 

The potential (|l|) 

(49) u^6p{x-fl; g2,g3) + 6p{x- '92,9^) 

with arbitrary invariants 92-, g^, 92, 93 is a finite-gap one for the equation (3), but 
its spectral characteristics can not be obtained in the framework of elliptic soliton 
theory. The corresponding commuting operator pencil is derived with the help of 
the equation (p2) and takes the form 



(9 (u - ci) A - 3 u'" + 6uu' + u') - 

- (27A2 + 9u'A -m('^> - 3u'2 + - ci u" - u'^ + 27c2) + 

-I- 6 A {u" - + u) * = /i 

Hence, the canonical form F{fi, A) = of the associated trigonal curve of genus 
9 — A is obtained by elimination of 

fj.^ +35(36a2 bA* - (a^ - 16g^ a* + 16 + 192 (b ~ 3 g.^) b + 192 ab'^) X"^ + 

-I- 48 (93 ~ 92 a'' b + 4 /X + 3^ {27 X"^ ~ 216 (a^b - 4 g.^ + 2 g^ b + 8b^) \^ - 

-2{a'^ + 30 92 a'^ -96g^a'^ -8 (45 fe^ _ 2I633 fe - 8g| + ^32 93) + 
+ 576 {b + 6g.^) g^o? b - 2^9 glab"^ - 2H^ (b^ -2g^ b^)) - 

- 288 (a^ - 2^2 + 24fe2)(g3 - 92 + 4b^) A) = 0, a = g^ - 92, b = g.^-g.^ 

and the corresponding ^'-function is given by the expression 

^, f XF^ -GH + FH' -F'H 

^) ^ '''''J g^-uF^-FH + F'G-FG' 
where prime denotes a derivation in x and 



F = ~3 u'" + 6 u u' - 3 ci u' + 9 (u - ci) A, ee 6 (w" - + u) A - ^, 

G = M<'^' + ci {u" + u2) _ 3 u'2 - 9 u' A - ^ u3 _ 27 (A2 + C2), 
93-93 8 (393 +93)92 - (393 + 93)92 



Ci = -12- . , _ „ 

92 - 92 ^ 3 92 -92 

One particular case of the potential ( p9[ ) and more general property of finite-gap 
potentials are discussed in |]39|] . 

The natural generalization of Hermite's method is to consider nonlinear homo- 
geneous ansatzs for the ^'-function. For instance, the quadratic ansatz 

^ ^ gfcx ^ $(a; _ Oj.; a) $(x - 0„; a). 

j, n 

However, this does not fit into the framework of finite-gap integration theory, be- 
cause the poles of the potential can depend on the spectral parameter. The following 
example with the transcendental dependence on spectral parameter elucidates this: 

*" = {6p{x)+2p{x~X)+'ip{X)) *(x;A) = ^^x; A). 
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Actually, the quadratic (and higher) ansatzs will not give an advantage due to the 
relation 

$(a;;a)2 = -$'(a;; 2 a, C(2 a) -2 ({a)) 
and we again arrive at the framework of Hermite's method. 

Note that the nonintegrability of equation (1) in context of the point 3) of the 
Proposition, nevertheless, can be useful for its integrability with special values 
of A or for more complex operator pencils with a dependence of the potential (say 
parameters a, b, c, d in Examples 4—5) on A. The availability of additional con- 
stants in the potentials may be considered as a family of spectral pencils, and under 
fixed values of A, as new spectral problems. For instance, the 2-gap Lame potential 
u — 6 pix) for equation (^ is obtained from example ( |3^ ) with A = and c? — s- A, 
VP' — > A less simple example is to swap the parameters A ^ c in equation (^6|), 
whereupon one finds the finite-gap operator A-pencil 

- 12 {p{x) + A^) - (c + 12 A p{x)) = 

with the algebraic curve F{k, A) = of genus g = 8 |, formula (15)]. 
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